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Plasma plumes with exotically segmented channel structure and plasma bullet propagation are produced in 
atmospheric plasma jets. This is achieved by tailoring interruptions of a continuous DC power supply over 
the time scales of lifetimes of residual electrons produced by the preceding discharge phase. These 
phenomena are explained by studying the plasma dynamics using nanosecond-precision imaging. One of 
the plumes is produced using 2 - 10 \is interruptions in the 8 kV DC voltage and features a still bright 
channel from which a propagating bullet detaches. A shorter interruption of 900 ns produces a plume with 
the additional long conducting dark channel between the jet nozzle and the bright area. The bullet size, 
formation dynamics, and propagation speed and distance can be effectively controlled. This may lead to 
micrometer- and nanosecond-precision delivery of quantized plasma bits, warranted for next-generation 
health, materials, and device technologies. 

Atmospheric-pressure non-equilibrium plasmas have recently received a major attention due to their 
widespread applications in surface and materials processing, biomedical engineering, and medicine 1 " 5 . 
Among them, plasma jets have been of particular interest owing to the highly-focused treatment 
spots, portability, near-room temperature of the plasma plume, as well as demonstrated efficacy and flex- 
ibility in a broad range of applications spanning from cancer treatment to nanotechnology 6 " 11 . The plasma 
jets feature a very unique dynamic behaviour compared to other types of atmospheric-pressure plasmas. 
Indeed, the discovery of fast-propagating plasma bullets (also termed pulsed atmospheric pressure streamers 
(PAPSs)) opens an exciting possibility to control small bits of plasmas over sub-microsecond and possibly 
even nanosecond time scales 12 " 16 . 

However, this fascinating phenomenon is still challenging to control and is far from being understood despite 
of a few mechanisms proposed 17 " 24 . Since these discharges are generated by time-varying electric fields, a series of 
breakdowns produces residual electrons within the channel of the plasma jet. The reported density of the residual 
electrons is vastly different by many orders of magnitude (e.g., from ~10 2 to ~10 2 cm" 3 ) 17,24 " 30 . Despite this 
quantitative uncertainty, it is generally believed that the residual electron dynamics plays a major role in the 
formation of the still plasma channel and the propagating plasma bullets that detach from it 19 . However, it is 
extremely difficult to control and monitor such a dynamic behaviour, mostly because of the very high (e.g., 
nanosecond) precision required as well as large uncertainties in the residual electron lifetimes. When using 
conventional photography, the bullet-prone plasma plumes commonly appear as dim glow channels. This is 
why despite the excellent recent progress in the studies of the plasma streamer phenomena, the understanding of 
the mechanisms of the effect of the residual electrons on the streamer propagation is far from being complete. 

To control the dynamics of the residual electrons, and hence, the plasma channel structure and the plasma 
bullet propagation, we customized the DC voltage to merely vary it over the timescales comparable with the 
expected lifetime of the residual electrons. In this way, by retaining the residual electrons long enough to affect the 
breakdown dynamics, here we report on the two exotic bullet-prone plume structures both featuring an uncom- 
mon large bright still area in addition to the common dim propagating bullet area. Moreover, one of these plumes 
also contains a highly- conducting dark channel behind the bright area. Importantly, these highly- unusual plasma 
plumes were generated using short interruptions of a continuous high-voltage DC supply, which is also a simpler 
and energy- efficient alternative to the commonly used short DC pulses. The observed effects are explained using 
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nanosecond time -resolved imaging of the discharge dynamics over 
the very short voltage rise and fall phases. This deterministic 
approach offers a simple and versatile control of the jet channel 
structure and conductivity as well as the plasma bullet propagation 
vitally needed for next-generation applications. 
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Figure 1 | Atmospheric plasma jet and DC pulse waveforms. 

(a) Schematic of the atmospheric plasma jet device, gas feed, power supply, 
and plasma diagnostics. Current-voltage pulsed DC waveforms 
corresponding to the pulse width (b) 2 (as, one single pulse is shown; 
(c) 998 (as, a 2 (as interruption between two pulses is shown; 
and (d) 999.1 (as, a 900 ns interruption between two pulses is shown. 
Breakdown occurs only during the voltage rise and fall phases. 



Results 

To demonstrate the highly-unusual plasma bullet-prone jet plumes 
by controlling the residual electrons, we used the atmospheric 
plasma jet device sketched in Fig. 1(a). The voltage waveforms used 
are shown in Figs. l(b)-(d). These waveforms were chosen to: (i) 
produce the breakdown merely during the rise and fall phases of the 
voltage pulse; and (ii) trigger the next voltage rise or fall phase with 
the delay comparable with the lifetime of the residual electrons pro- 
duced during the previous discharge phase. Using a typical observ- 
able plasma bullet propagation time of 1 (as as an (assumed) estimate 
of the residual electrons lifetime, the duration of such a delay was 
varied from 900 ns to 2 (as. 

Fig. 1(b) shows a single 8 kV DC voltage pulse of a 2 (as duration 
(and 1 ms repetition rate) and the current it produces during the 
initial rise phase and the following fall phase. The continuous 8 kV 
DC voltage was also repeatedly interrupted for 2 (as, also every 1 ms, 
as shown in Fig. 1(c). The interruption time was also reduced to 
900 ns retaining the same repetition rate (Fig. 1(d)). During the 
interrupted DC operation, the fall phase preceded the rise phase as 
shown in Fig. l(c,d). In all cases, the discharges were only produced 
during the rise or fall phases as evidenced by the current spikes in 
Fig. l(b-d). This was also chosen deliberately to clearly separate the 
different discharge phases during the high-speed, time-resolved 
imaging and focus on the plasma bullet propagation which was prev- 
iously reported to merely occur during the rise phase 13 . 

Remarkably, all the three waveforms in Fig. l(b-d) produced three 
totally different plasma plumes depicted in Fig. 2(a-c). The 2 (as 
pulse has generated a common plume with a large dim area indicative 
of a propagating bullet (Fig. 2(a)). Importantly, this type of the 
plasma plumes was generated using DC pulses with the duration 
up to 990 (as at the same repetition rate of 1 ms. However, the other 
two interrupted DC waveforms produced exotic plumes with two 
and three distinctive areas labelled in Fig. 2(b) and (c), respectively. 
When the break between the pulses was 2 (as (Fig. 2(b)), the plume 
appears as a bright channel (area I) and a dim (likely bullet propaga- 
tion) area II. This type of structures was revealed for the 2 - 10 (as 
interruptions of 8 kV DC pulses. 

When the duration of DC interruption was reduced to less than 
1 us, a very different plume structure was discovered (Fig. 2(c)). One 
can clearly see the dark area I, the bright area II, and the dim (likely 
bullet propagation) area III in the plume generated using a 900 ns 
interruption (Fig. 1(d)). 

To help explaining the observed exotic plume structures, time- 
resolved, high-speed imaging was implemented selectively during 
the rise and fall phases of each of the waveforms in Fig. 1. In all cases, 
it was found that the plasma bullets are only produced during the rise 




Figure 2 | Photographs of three distinctive plasma plume types generated 
using, (a) 2 (as pulse; mostly a dim (typical to a plasma bullet) area is seen; 

(b) 2 (as interruption; bright (I) and dim (II) areas are seen; and (c) 900 ns 
interruption; dark (I), bright (II), and dim (III) areas are seen. Images (a)- 

(c) produced using a 1 (as exposure from the beginning of the voltage rise 
phases labelled "0" in Figs. l(b)-l(d), respectively. 
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phases, for both the short pulses and the short interruptions. These 
high-speed images for the 2 jas pulse and the 2 |is and 900 ns breaks 
are shown in Figs. 3-5, respectively. The beginning of the time count 
was synchronized with the onset of the voltage rise phase. From these 
images one can clearly spot the very different discharge dynamics in 
each case. 

When a short 2 |is pulse is used, the plume dynamics (Fig. 3) is 
typical to the plasma bullet propagation reported previously 13 . The 
bright plasma channel originates within the nozzle and starts exiting 
it at 528 ns into the pulse. At 628 ns, the plume elongates to ~ 1 cm 
forming a bright bullet at its head. 100 ns later, the bullet separates 
from the rest of the plume. Other images in Fig. 3 show that the bullet 
propagates further, with the speed of the ionization front ~10 5 m/s 
(note that the speed of the gas flow is much lower, typically of the 
order of —10 m/s), eventually slowing down and shrinking in size. 
The rest of the plasma plume darkens significantly while the bullet 
propagates. 

Fig. 4 shows a very different dynamics of the plasma plume gen- 
erated by a short 2 |is interruption. The plume brightens and exits 
the nozzle much faster thereby forming a bright, non -propagating 
channel at 278 ns which then fully exits from the nozzle 100 ns later. 
The length of this channel is approximately 6 mm. At 428 ns, the 
channel brightens and the plasma bullet forms and tends to detach 
from the channel. 50 ns later, the bullet is fully detached and a dark 
channel forms between the head of the bright plasma channel and the 
small bullet-like plasma. As the bullet propagates further, its speed 
and size decrease while the luminosity of the bright channel I in 
Fig. 2(b) decreases. The speed and the length of the bullet propaga- 
tion are noticeably smaller than in Fig. 3 thereby leading to a shorter 
plasma plume in Fig. 2(b) compared to Fig. 2(a). For example, the 
propagation speed of the bullet is approximately 3 times smaller than 
in the 2 jas pulse case. 
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Figure 4 | Same as in Fig. 3 for the 2 \is DC interruption. The time labeled 
on each photograph corresponds to the time in Fig. 1(c). 

In order to elucidate whether the bright plasma channel I in 
Fig. 2(b) is formed locally or is due to the integrated effect of a fast 
moving bullet, a dedicated experiment using the ICCD camera with 
an exposure time of 2 ns and a step of 2 ns was performed. This high- 
speed imaging confirms that the bright plasma channel at 378 ns is 
clearly not due to the integrated effect of a fast moving bullet. 

If the break between the voltage fall and rise phases is further 
reduced, the bright plasma channel starts forming within the nozzle 
and then moves forward much faster. This can be seen in Fig. 5 which 
corresponds to the waveform in Fig. 1(d). Already at 173 ns, the 
bright channel clearly separates from the nozzle and leaves an 










728ns 




Figure 3 | High-speed photographs of the plasma plume generated 
during the voltage rise phase of the 2 ^s pulse. The exposure time is 5 ns. 
The time labeled on each photograph corresponds to the time in Fig. 1(b). 




Figure 5 | Same as in Fig. 3 for the 900 ns DC interruption. The time 
labeled on each photograph corresponds to the time in Fig. 1(d). 
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approx. 4 mm dark channel behind. The length of the dark channel 
between the nozzle and the left end of the bright plasma channel 
becomes smaller when the voltage interruption time decreases. At 
about 273 ns, a plasma bullet is formed at the right end of the bright 
plasma channel and starts propagating forward. 100 ns later, a clear 
gap between the bullet and the channel is visible. At this moment, the 
channel becomes dimmer and this tendency continues thereafter. 
Interestingly, the length of the dark space (area I in Fig. 2(c)) between 
the nozzle and the left end of the bright channel does not noticeably 
change with time. As the bullet propagates, it also reduces in size and 
slows down, similar to the two other cases. The propagation speed of 
the bullet in Fig. 5 is somewhat smaller than in Fig. 4. 

Discussion 

Here we discuss the observed phenomena and relate them with the 
density and response of residual electrons, in the context of the 
existing knowledge on the propagation mechanism of pulsed atmo- 
spheric plasma streamers. We will then summarize the interpreta- 
tions of the key observations and the remaining open issues. 

Dark channel. Numerical modelling shows that electron density of 
dark channel left behind the plasma bullet is typically of the order of 
10 11 - 10 12 cm" 3 17 . The high electron density results in a high 
conductivity of the dark channel. Thus, the electric field within the 
dark channel is low, typically about a few kV/cm 17 . Hence, the dark 
channel effectively serves as an extended electrode. This is also the 
reason why the plasma bullet propagation is inhibited at the falling 
edge of the applied voltage 18 . Because the conductivity of the dark 
channel is high, the plasma bullet is electrically connected to the 
high-voltage electrode. 

The electron density in the dark channel in Fig. 3 decreases with 
time and distance from the nozzle. The electron density in the dark 
channel behind the plasma bullet in Fig. 5 is also expected to behave 
quite similarly because of the similar trends of the electron density 
decay, as in the case of the pulse width of 2 |is. Quantitatively, the 
electron densities in the two dark channels in Fig. 5 are quite different 
over the most of the pulse cycle. The dark channel behind the bright 
channel has a lower mole fraction of air than the dark channel behind 
the plasma bullet. Consequently, the electron density behind the 
bright channel experiences a slower decay than behind the plasma 
bullet. More importantly, the electrons in these two dark channels are 
generated at different time. While the electrons in the dark channel 
behind the bright channel are generated during the discharge at both 
the falling and the rising edges of the pulse, the electrons in the dark 
channel behind the plasma bullet are generated during the discharge 
at the rising edge of the pulse (the discharge at the falling edge does 
not reach this position as can be seen in Fig. S3). It is estimated that 
the electron density before the falling edge discharge in the dark 
channel behind the bright channel is of the order of 10 9 cm" 3 , while 
in the dark channel behind the plasma bullet it is of the order of 
10 4 cm" 3 . This difference is due to the different 0 2 concentration. 
After the falling edge discharge, the electron density in the dark 
channel behind the bright channel rise to ~10 12 cm" 3 , while the 
electron density in the dark channel behind the plasma bullet 
remains much lower (of the order of 10 4 cm" 3 ) because the plasma 
did not reach the position. For the pulse width of 999.1 us, the 
discharge at the rising edge is ignited only several hundred ns after 
the discharge at the previous falling edge. Then the electron density 
in the both dark channels could reach the value of the order of 
~10 12 cm" 3 . 

2-10 \is breaks between voltage phases. In our experiments, when 
the voltage interruptions are longer than 2 /us but shorter than 10 us, 
a bright channel is generated in front of the nozzle at the rising edge. 
A plasma bullet is generated at the right tip of the bright channel and 
moves forward. This phenomenon does not appear when the voltage 



break is longer than 10 |is. It is thus reasonable to assume that the 
observed unusual plasma plume structure is due to the increase of the 
residual electron density left from the discharge at the falling edge of 
the previous voltage pulse. When the time gap between the rise phase 
and the preceding fall phase becomes snorter so does the break 
between the two consecutive discharges ignited during the both 
phases. In this case, species produced by the fall-phase discharge 
may not completely recombine before the onset of the subsequent 
discharge at the voltage rise phase. Therefore, because of the short 
time between the two discharges, the density of residual electrons 
produced during the preceding fall phase may remain high before the 
next discharge is ignited at the rising edge of the voltage pulse. We 
emphasize that the length of the plasma plume during the fall phase is 
much shorter compared to the rise phase, as can be clearly seen from 
Fig. 4 and Fig. S3(b). From Fig. 2(b) one can also note that the length 
of the bright section I is close to the length of the plasma channel 
deduced from high-speed imaging in Supplementary Fig. S3. It is 
thus reasonable to assume that the residual electron density within 
section I is relatively high. 

Consequently, section I becomes conductive at the onset of the rise 
phase. The dark channel is left from the propagation of the plasma 
bullet. The reason why section I appears bright is probably because 
the electric field within section I is higher than within the dark 
channel. The relatively high electric field within section I sustains a 
certain degree of ionization and excitation, which is balanced by the 
recombination, attachment and diffusion of the gas species. This is 
why the plasma of section I remains stable for a relatively long time as 
evidenced in Fig. 4. For section II shown in Fig. 2(b), the residual 
electron density is much lower than in section I which helps sustain- 
ing the plasma bullet propagation as seen in Fig. 4. The numerical 
modelling of the distribution of the residual electron density is pre- 
sented in Sec. S4 of Supplementary Information. As shown in Fig. S7, 
when the pulse width is 998 |is, the residual electron density is about 
4 X 10 11 cm" 3 in section I in Fig. 2(b) and about 10 4 cm" 3 in section 
II in Fig. 2(b). This result is consistent with the phenomenological 
discussions in this work. 

Shorter (<2 \is) breaks between voltage phases. When the voltage 
interruption is shorter than 2 |is, the residual electron density in 
front of the nozzle is even higher when the discharge is ignited 
during the rise phase. Hence, section I in Fig. 2(c) has an even 
higher electric conductivity and a weaker electric field. The weaker 
electric field results in a low electron temperature, so the electrons 
cannot gain enough energy to ionize or excite atoms or molecules. 

Role of air diffusion. It should be pointed out that the diffusion of air 
plays a significant role in the phenomena discussed above. Due to the 
diffusion, the air concentration becomes higher and higher in the He 
stream with the increase of the distance from the nozzle. Because of 
the very high electron affinity of oxygen, the high concentration of air 
in the He flow results in the fast decay of the residual electron density. 
This is why the residual electron density decreases away from the 
nozzle. Consequently, sections I and II in Fig. 2(c) feature the high 
and medium residual electron density. This in turn results in the dark 
region I, and bright region II, respectively. As shown in Figure S7, 
when the pulse width is 999.1 jas, the residual electron density 
decreases with the distance from the nozzle in sections I and II in 
Fig. 2(c) (see Sec. S4 in Supplementary Information). This numerical 
result further supports our discussion. 

Moreover, diffusion of oxygen and nitrogen species contributes to 
de-excitation of He metastables, which contributes to the optical 
emission, see the optical emission spectra in Figure S5 in Supple- 
mentary Information. The mechanism is related to the Penning ion- 
ization whereby these molecules are ionized. It is well known that 
Penning ionization can enhance the discharge. It increases the elec- 
tron density and the velocity of the plasma bullet. This is why this 
effect maybe considered as counteracting to the effect of the decrease 
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of the residual electron density. The resulting optical emission from 
the plasma plume is a consequence of the ionization-excitation-radi- 
ation balance. The residual electron- and helium metastable-related 
mechanisms both contribute to this balance, and consequently, to the 
observed discharge behaviour. 

Summary of observations and interpretations. We have thus 
observed two new phenomena produced when the duration of 
interruption of DC voltage are comparable with the lifetime of the 
residual electron density. First, when the break is between 2 |is and 
10 |is, the plasma jet generated during the voltage rise phase can be 
divided into two sections. In section I which appears next to the 
nozzle, the plasma appears like a bright channel and has a length 
of several millimetres. The bright channel appears in section I and a 
plasma bullet is generated at the tip of the bright channel. The bullet 
then detaches and propagates along the He stream. The bright 
channel is due to the high residual electron density left from the 
discharge at the preceding voltage fall phase. The high residual 
electron density results in the relatively high electric conductivity 
and hence a weak electric field. The electron density of the bright 
channel of section I is lower compared to the dark channel due to 
recombination, attachment, and other loss effects. Consequently, the 
electric field within section I is stronger than in the dark channel. 
This electric field leads to the higher electron temperature which in 
turn raises the level of excitation of atomic/molecular species in the 
gas flow. This produces intense optical emission from section I. For 
section II in Fig. 2(b), the residual electron density is much lower 
compared to section I and the plasma retains the propagating bullet 
behaviour. 

Second, when the voltage break is further reduced, the gap 
between the discharge at the rise phase and the discharge at the 
preceding fall phase becomes shorter and shorter, which increases 
the residual electron density in the channel. The bright plasma chan- 
nel (section II in Fig. 2(c)) detaches further and further away from the 
nozzle as the voltage break decreases. This creates a dark channel 
(section I in Fig. 2(c)) between the nozzle and the left end of the 
bright plasma channel. On the right end of the bright plasma chan- 
nel, the jet still retains the plasma bullet behaviour. Due to the air 
diffusion into the He stream and the discharge at the falling edge of 
the previous voltage pulse only reaching sections I and II, the residual 
electron density in sections I, II, and III of Fig. 2(c) is high, medium, 
and low, respectively and the residual electron density in sections I 
and II are many orders higher than that in section III. This difference 
in the electron densities results in the appearance of the dark, bright, 
and bullet areas, respectively. The numerical modelling of the distri- 
bution of the residual electron density presented in Sec. S4 of 
Supplementary Information supports this explanation. 

Therefore, our results suggest that the residual electron density in 
the discharge channels of atmospheric-pressure plasma jets plays a 
critical role in the formation, structure, and temporal dynamics of the 
plasma plumes. This phenomenon has not been studied previously 
despite a large number of studies on the effect of seed electrons on the 
streamer propagation. Precise tailoring of the duration and ramping 
of DC power supply appears to be a viable way to achieve the exotic 
multi- segmented plumes and control the plasma bullet formation 
and propagation. Interaction of the plasma jet with the surrounding 
air is also a decisive factor in the spatial distribution and dynamics of 
the residual electron density. Detailed study of this interaction and 
other possible mechanisms is therefore highly warranted in the near 
future. For example, the absorption of photoemission by He, N 2 and 
0 2 may also play a significant role in the streamer propagation 31 . 

Potential applications. The demonstrated control of the plasma 
bullet formation and propagation is highly-promising for the 
controlled delivery of quantized bits of plasmas with micrometer 
and nanosecond precision, which may open new exciting prospects 
for applications in cellular-level medical therapies and selective 



treatment of microscopic objects in numerous materials and device 
technologies. Indeed, although the typical size of the spot of the 
plasma jet is in the sub -mm range, it can be reduced to a couple of 
hundreds or even a few tens of micrometers by a proper optimization 
of the discharge parameters. Apart from the small dimensions of 
the localized treatment spots, the surface temperature within these 
spots can reasonably be expected to be close to room temperature. 
This is particularly important for the synthesis and processing of 
temperature-sensitive materials such as polymers 32,33 . Moreover, 
low energies of the plasma ions (the ion velocities in front of the 
plasma bullet are of the order of ~ 10 3 m/s, which is much lower than 
typically in low-pressure plasmas) are very favourable for the 
treatment of materials that are highly- susceptible to ion damage 
such as delicate graphenes and single-walled carbon nanotubes 34 ' 35 . 
These interesting features of the pulsed atmospheric pressure 
streamers may bring several advantages to the field of nanoscale 
synthesis and processing where low-pressure plasmas have demon- 
strated a number of advantages over many other conventional 
methods 36 " 38 . 

Methods 

Fig. 1 (a) is the schematic of the experimental setup. The plasma plume is generated by 
a single electrode plasma jet device. The voltage pulse repetition rate (ON + OFF 
pulse sequences) is fixed at 1 kHz throughout this paper. The pulse width was varied 
from 900 ns to 999. 1 us. The applied voltage is measured by a P6015A Tektronix HV 
probe and the current by a Pearson 6585 current probe. A fast ICCD camera 
(Princeton Instruments, Model: PIMAX2) is used to capture the nanosecond 
dynamics of the discharge. The exposure time is fixed at 5 ns for the photographs 
shown in this paper. A 1 us exposure is used to produce still images of the plasma 
plumes in Fig. 2. The terms "dim", "dark", and "bright" are used to describe the 
relative luminance in different areas of the plasma jet rather than to express their 
absolute luminance. Similar presentations of the relative luminance such as "dark 
space", "glow" are consistently used to describe the visual appearance of the glow 
discharge. The discharge currents shown in Fig. 1 are the total currents and represent 
the sum of the displacement current and the actual discharge current, where the 
actual discharge current is much smaller than the displacement current. It should be 
mentioned that the gas temperature may be affected by the pulse repetition frequency; 
however, this frequency was fixed at 1 kHz in our experiments. We have also 
observed that when varying the pulse width, the discharge current also changes; 
however, this change is fairly small. As a result, the gas temperature does not undergo 
any significant changes. During the jet operation at a fixed pulse repetition frequency, 
the gas temperature remains close to room temperature. This is why the assumption 
of the constant near-room temperature is reasonable. For each waveform in Fig. 1, the 
zero time moment is set at the beginning of the rise phase. The reason for this choice 
was to temporally resolve the plasma bullet propagation which was only observed 
during the voltage rise phase. The gas flow used in this series of experiments was 1 
1/min. Under such conditions, the gas flow is laminar and the pulsed atmospheric 
pressure streamer propagates in a linear fashion. Above ~5 1/min, turbulent gas flows 
develop, which make the streamer propagation more sophisticated. For example, it 
was observed that the plasma bullet may exhibit a snake-like propagation, which was 
intentionally avoided for the purpose of this work. More details about the plasma jet 
device, the associated measurements and the method used to calculate the electron 
density can be found in Supplementary Information. 
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